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The effects of turbulent intensity and vortex scale of simulated natural wind on performance of a horizontal axis 


wind turbine (HAWT) are mainly investigated in this paper. In particular, the unsteadiness and turbulence of wind 


in Japan are stronger than ones in Europe and North America in general. Hence, Japanese engineers should take 


account of the velocity unsteadiness of natural wind at installed open-air location to design a higher performance 


wind turbine. Using the originally designed five wind turbines on the basis of NACA and MEL blades, the de- 


pendencies of the wind frequency and vortex scale of the simulated natural wind are presented. As the results, the 


power coefficient of the newly designed MEL3-type rotor in the simulated natural wind is 130% larger than one 


in steady wind. 
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Introduction 


A wind turbine is developed through the following 
procedure in general. Firstly, a design developer eva- 
luates performance of a newly designed wind turbine in 
steady wind through a wind tunnel test to determine the 
basic design. Secondly, performance of a scaled model of 
a prototype wind turbine is investigated at the installation 
location for a year to evaluate the performance. Finally, 
the developer determines the design. For these proce- 
dures, it needs long time and large cost to develop a new 
wind turbine. 


Received: October 2016 


Recently, it is important that one should consider the 
effects of unsteadiness and turbulence of the wind on the 
performance for high output power. This is because that a 
wind turbine in natural wind with velocity fluctuation 
and turbulence show larger power than one in steady 
wind for some installation locations [1]. For the reasons, 
the current wind turbine design is developed in consider- 
ation of disturbance effect of the wind [2, 3]. 

The purpose of this study is to establish the method of 
a wind turbine design, which is tailor-made for the natu- 
ral wind at installation location without the field test. In 
order to evaluate the performance in natural wind, simu- 
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Nomenclature 
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U,u Wind velocity (m/s) 
Greek letters 
P Fluid density (kg/m?) 
oO Standard deviation of wind velocity 
w Circular frequency (rad/s) 
a Tip speed ratio 
Subscripts 
D Design value 
S Steady wind, steady component 
u Unsteady wind, unsteady component 
© Time averaged value 
O Unsteady value 
Torque Detector DC Motor 


(Maxon RE25) 


Steady Flow } Rotor~, : 
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Fig. 1 Schematics of experimental system to measure wind turbine performance. 


Table 1 Design condition of rotor 


Kott Up a AOA POT MOT Mass 
[m/s] 2 [deg.] (x./c) [kgm'] [g] 

NACAI 10 5.0 - 50% 1.52x10* 120 
NACA2 6.5 3.5 8 35% 8.65x10° 102 
MEL1 6 5.0 10 33% 2.81107 67 
MEL2 6 4.0 15 33% 3.95x 10° 71 
MEL3 6 2.5 15 33% 1.42x10* 104 


Ap=Designed Tip speed ratio, AOA=Angle of attack, POT=Position of torsional axis from leading edge, MOI=Moment of inertia. 


lated natural wind is generated by a multi-fan type active 
control wind tunnel. In this study, output powers of five 
rotors are examined in the simulated natural wind which 
is made on the basis of Karman’s power spectrum density 
function. The effect of vortex scale on the wind turbine 
performance and appropriate design point are shown [4]. 


Experimental Procedure 


Experimental apparatus 


Figure 1 shows the measuring system in experiment. 
The wind turbine type is a horizontal axis with three 


blades. The cross section of the wind tunnel is 1mx1m. 
The flow is generated by the actively controlled 66 fans. 
The wind turbine power measuring system is consisted of 
the torque detector, the rotational speed sensor, the torque 
converter (Onosokki, SS-005, MP-981), the DC motor 
(Maxon, RE25, 20W) and the DC electric power supply. 
The DC motor works as a power generator to control the 
rotor load and the rotational speed. 


Wind Turbine Models 


Table 1 and Figure 2 show the design conditions of 
rotors, its schematic and photo respectively. The rotors of 
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NACA are designed on the basis of NACA 632xx series 
blade of Reynold number Re = 3.0 x 10°. The rotor NACA1 
is designed by Furukawa et al. for a long type wind lens 
turbine [5], which the profile is formed from NACA63218 
at root to NACA63212 at tip along a span. The characte- 
ristics and performance for steady and sinusoidally oscil- 
lating velocity wind have been studied in the authors’ 
previous papers [6]. The rotor NACA2 and MEL1-3 are 
newly designed on the basis of NACA63218 with Re = 
3.0 x 10° and MEL002 with Re = 9.0 x 10* respectively. 
Here, the MEL series blade (MEL002) was developed by 
the National Institute of Advanced Industrial Science and 
Technology (AIST) for wind turbine in Japan. 


Fig. 2 Schematic and photograph of the wind turbine model, 
MEL3. 


In this experiment, two kinds of Reynolds number 
should be considered. One is Rotor Reynolds Number, 
Re=9.0x10*, which is calculated by steady velocity 7 m/s 
and rotor diameter 0.194 m as a characteristic length. 
Another is local blade Reynolds number. For example, 
when inflow velocity of the blade is assumed as 2 m/s, it 
is determined as Re = 1.8 x 10*~ 3.7 x 10* with tip speed 
ratio 2 = 2 ~ 4 and blade chord length 2cm as a characte- 
ristic length. Then, the local blade Reynolds number is 
the same order as one of MEL002. 


Definition of Power and Torque Coefficients 


The power and torque coefficients of wind turbines are 
defined by following equations; 


G= (1) 
— pAU? 
2 
(oe (2) 
— pArU” 
ha 
A= (3) 


Here, P,T,p,Aand U denote generated power and 
torque of a wind turbine, fluid density, rotor swept area 
and the upstream time averaged flow velocity, respec- 
tively. 


Wind Turbine Performance in Steady Flow 


Figure 3 and 4 show the effect of rotor type on the 
power and torque coefficients in steady wind of 7 m/s. 


The rotors NACA1 and MEL3 show larger power coeffi- 
cients 0.32 than other rotors. The maximum power coef- 
ficient of NACA2 and MEL 1-3 exist at à = 3.0, 3.8, 4.1 
and 2.3, respectively. The tip speed ratios are more of the 
same as the design tip speed ratios. Thus, it means that 
the authors’ design is appropriate in this study. 
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Fig. 3 Power coefficient of the wind turbines in steady wind 
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Fig. 4 Torque coefficient of the wind turbines in steady wind 
at U =7m/s 


Performance in Sinusoidally Oscillating Velocity 
Flow 


Wind conditions 

The upstream wind velocity of the wind turbine is un- 
dergoing identical harmonic oscillated as follow: 

U(t)=U +i sin(2zft) (4) 

The performances of the wind turbines are investi- 
gated as the upstream mean velocities U =7m/s with 
oscillating amplitudes wu =1.4m/s (20% of 7m/s). The 
wind velocity oscillates at frequencies f = 0.033, 0.05, 
0.083 and 0.25Hz and their periods T,,= 30s, 20s, 12s, 8s, 
4s. Here, the amplitude, frequencies and periods are de- 
termined on the basis of the natural wind data of Kasuga 
campus of Kyushu University in Japan [7]. 


Increasing ratio of power coefficient 


The maximum power coefficient ratio Cy,/Cy; is de- 
fined to compare the performances with unsteady and 
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steady winds. Here, C,, and Cy, are maximum power 
coefficients in the steady and the sinusoidally oscillating 
velocity wind respectively. Figures 5 and 6 show the ratio 
Cwu/Cws at A = 2.0 (1400 min”) and 2 = 3.9 (2700 min”). 
MEL3 has the large increasing ratio for all period Tw at A 
= 2.0 (1400 min’) and T,<12s at A = 3.9 (2700 min’). In 
particular, its largest increasing ratio is 130% at Tọ =20s 
and A = 3.9. While C,,,/C,, of NACAI1 and MELI! are 
larger than one of MEL3 at T, =4s and å = 3.9 (2700 
min`’). As the results, one can conclude that MEL3 is 
most appropriate rotor for the oscillating velocity wind. 
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Fig. 5 Power coefficient ratio with wind turbine type and 
wind period in the oscillating velocity wind at \=2.0. 
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Fig. 6 Power coefficient ratio with wind turbine type and 


wind period in the oscillating velocity wind at \=3.9. 


Performance in Simulated Natural Fluctuating 
Velocity Flow 


The simulated natural wind with fluctuating velocity 
based on Karman’s power spectrum density function 
(PSD) of Eqs. (5) and (6). 


S,(f) =4V7LU 


376 (5) 


(6) 
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Here, I and L mean turbulent intensity and vortex scale 
respectively. The velocity data of time history are gener- 
ated through inverse Fourier transformation with random 
phase in this study. 


Wind conditions 


Table 2 shows the conditions of the simulated natural 
wind to examine the effect of I and L on power coeffi- 
cient. The experimental duration time is 10 minutes for 
each case. The mean velocity, I and L of the experimental 
simulated wind do not accurately correspond to target 
conditions. Hence they are monitored during measure- 
ment by hot-wire anemometer. The data show the column 
of “Measurement” in Table 2. 

Figure 7 presents the comparison between theoretical 
data and experimental result of PSD [S,(f)] at target con- 
ditions U = 6.5m/s, I = 5.0% and L = 6 m. The theoreti- 
cal and experimental data is almost the same for the fre- 
quency 0.03 < f < 1.0. It means that the simulated wind is 
appropriately generated on the basis of Karman’s PSD of 


Eq. (5). 


Table 2 Conditions of turbulent wind. 


Target Measurement 
U I L U I 
[m/s] [%] [m] [m/s] [%] 
3 6.8 4.5 
NACA1 6.5 5 6 6.4 5:1 
10 7.1 5.5 
3 6.5 4.5 
MEL3 6.5 5 6 6.4 5.0 
10 6.3 5.4 
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Fig. 7 Comparison of power spectrum between theory and 
experiment for fluctuating wind velocity. 


Experimental results and discussions 


Figure 8 and 9 show the effect of vortex scale L on the 
power coefficients of NACAI and MEL3 with the target 
turbulent intensity I = 5% in Table 2. In addition, power 
coefficient in steady wind of 7m/s is also presented in the 
figures as reference. The gray region of “resonance” 
show that rotor axis mechanically resonates with rotation. 
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Hence the data of resonance region are omitted. 

The maximum power coefficient of MEL3 is the same 
as one of NACA1. The maximum power coefficient ra- 
tios (Cy,/Cys) of MEL3 and NACAI are 1.30 and 1.16 
respectively. 

The power coefficient of NACA1 at L = 6 m is almost 
the same as one at L = 10 m. Thus, the effect of L is not 
large on the power coefficient of NACAI at I = 5%. 
While the power coefficient of MEL3 increases with in- 
creasing L. Hence, the effect L is large on the power 
coefficient of MEL3 at I = 5%. The maximum C,,, is oc- 
curred with L =10 m. This is because that local blade 
Reynolds number is the same order as Reynolds number 
of MEL-type blade. Thus, it seems that MEL3 rotor is 
responsible for the vortex scale. According to the results, 
MEL3 shows the highest performance for the natural 
winds. 
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Fig. 8 Effect of vortex scale on power coefficient of NACA1 
in the condition in Table 2. 
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Fig.9 Effect of vortex scale on power coefficient of MEL3 in 
the condition in Table 2. 


Conclusions 


The paper presents the performances of the five pro- 
totype wind turbines for the unsteady fluctuating velocity 
wind — sinusoidally oscillating velocity wind and simu- 
lated natural wind. The rotors with four blades are newly 
designed on the basis of NACA five number series 


(NACA63xxx) and MEL002 blades. 

According to the experimental results, the conclusions 
are as follows. 

1. It is better that the designed tip speed ratio is in lower 
range in the sinusoidally oscillating velocity wind. 

2. The rotor MEL3 shows the highest performance in the 
five rotors for three kind of winds - steady, sinusoi- 
dally oscillating and simulated natural winds. In par- 
ticular, the maximum power coefficient in the simu- 
lated natural wind is as 130% large as one in the 
steady wind. 

3. Wind turbine performance can be increased through 
tailoring wind rotor design to turbulent character of 
installing location wind. 
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